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ABSTRACT 

Highly ionized, z = 0 metal absorption lines detected in the X-ray spectra of 
background active galactic nuclei (AGNs) provide an effective method to probe the hot 
(T r\j 10 6 K) gas and its metal content in and around the Milky Way. We present an 
all-sky survey of the K a transition of the local O VII absorption lines obtained by Voigt- 
prohle fitting archival XMM-Newton observations. A total of 43 AGNs were selected, 
among which 12 are BL Lac-type AGNs, and the rest are Seyfert 1 galaxies. At above 
the 3cr level the local O VII absorption lines were detected in 21 AGNs, among which 
7 were newly discovered in this work. The sky covering fraction, defined as the ratio 
between the number of detections and the sample size, increases from at about 40% for 
all targets to 100% for the brightest targets, suggesting a uniform distribution of the 
O VII absorbers. We correlate the line equivalent width with the Galactic coordinates 
and do not find any strong correlations between these quantities. Some AGNs have 
warm absorbers that may complicate the analysis of the local X-ray absorber since the 
recession velocity can be compensated by the outflow velocity, especially for the nearby 
targets. We discuss the potential impact of the warm absorbers on our analysis. A com¬ 
prehensive theoretical modeling of the X-ray absorbers will be presented in a later paper. 


Subject headings: quasars: absorption lines 


1. Introduction 


Since the launch of the Chandra and XMM-Newton X-ray Observatories, a number of z ~ 0 
absorption lines, produced by highly ionized metal species, were detected in the spectra of back¬ 


ground active galactic nuclei (see, e.g., Nicastro et al. 2002 Fang et al. 2002 Rasmussen et al 

at 

~ 2 imply a hot, metal enriched gas distribution in and around 


2003). These highly ionized metal absorption lines, mostly in the form of O VII and with a column 
density on the order of ~ 10 16r 


cm 
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ALL SKY MAP 



Fig. 1.— All-sky Hammer-Aitoff projection of our targets in the Galactic coordinates. Red squares 
are BL Lac-type targets, and blue circles are Seyfert 1 galaxies. 


our Galaxy. It stirs a great amount of interest for the potential impact on our understanding of 
galaxy formation and evolution. 


However, ever since the first detection, a crucial debate is the origin of the X-ray absorbers. For 
the local absorption lines seen in the quasar spectra, even with the high resolution spectrometers 
on-board Chandra and XMM-Newton, we can only limit the location of the absorbers to within 
a few Mpc, therefore it is still unclear whether this gas is associated with our disk, an extended 
galactic halo (Wang et al. 2006; Fang et al. 2006), or even permeates the Local Group in the 
form of the intragroup medium (e.g., Nicastro et ah 2002; Williams et ah 2005). Identifying 
the origin of this gas has significant implications on many fronts. A disk origin would benefit the 
study of the three-phase model of the interstellar medium (ISM, Ostriker & McKee 1978). If the 
absorption is produced by the hot, circumgalactic medium (CGM) in the distant halo, it would 
imply a significant reservoir of baryons may account for the “missing Galactic baryons” (see, e.g., 


Fang et ah 2006; Gupta et ah 2012). Furthermore, the hot gas could reside in the hot intergroup 


medium in the Local Group, which would potentially explain the so-called warm-hot intergalactic 
medium (WHIM, see, e.g, Nicastro et ah 2002). 


The existence of the hot gas in and around the Milky Way has long been established through 
the observations of the soft X-ray background emission (see, e.g., McCommon Sz Sanders 1990; 
Wang Sz Yu 1995; Snowden et ah 1998; Kuntz & Snowden 1998; Galeazzi et ah 2007; Henley Sz 
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Fig. 2.— Redshift distribution of the sample AGNs. 


Shelton 2008, 2013; Yoshino et al. 2009). However, current data cannot distinguish between a disk 
and a halo origin of the X-emitting gas. While most studies have argued for a disk morphology of 


the hot gas (see, e.g., 

Yao & Wang 

2007, 

Hagihara et al. 

2010 

), other studies have suggested that 

a halo origin is also possible ( 

Fang et al. 2013 

Henley & Shelton 2013 

)■ 


Substantial progress has been made since these absorbers were first identified. McKernan 


et al. (2004) studied the Chandra spectra of 15 type I AGNs and argued that the observed metal 


absorption lines should be local to our Galaxy instead of high velocity outflows intrinsic to the 
AGNs. Fang et al. (2006) further investigated 20 AGNs from both Chandra and XMM-Newton 
observations and concluded the absorption line systems must be in/around the Milky Way instead 


of in the intragroup medium of the distant Local Group. Bregman & Lloyd-Davies (2007) reached 


a similar conclusion by studying the correlation between the X-ray absorption line strength and the 
ROSAT background emission measure. Furthermore, based on a survey of 26 AGNs, Bregman & 


Lloyd-Davies (2007) and Miller & Bregman (2013) analyzed the structure of the CGM and argued 


that the hot gas accounts for at most 10-50% of the missing Galactic baryons. 

Despite all the progress, a critical, unresolved issue is still the relative contribution of the disk 
and the distant halo to the observed absorption seen in X-ray. Inspired by the disk morphology 
of the X-ray emission of nearby galaxies, several authors (see, e.g., Yao et al. 2009; Hagihara et 
al. 2010) argued that most of the X-ray absorbing gas is confined within a few kpc of the Milky 
Way disk. X-ray absorption lines were also detected in the nearby Galactic X-ray binaries (see, 


e.g., Miller et al. 2004; Yao & Wang 2005). Such absorption lines have been compared with those 


detected in the AGNs to determine the relative contribution of the disk and halo (see, e.g., Wang 
et al. 2005; Yao & Wang 2007). However, A complication arises in these studies because the X-ray 
absorpriton lines in the XRBs may be contaminated by the circumstellar medium intrinsic to the 


XRBs (Miller et al. 2004 Cackett et al. 2008). While a disk origin for the hot gas is attractive, a 
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halo-origin cannot be definitely excluded (see, e.g., Fang et al. 2006; Gupta et al. 2012; Miller & 


Bregman 2013, but also see Wang Sz Yao 2012) 


Analysis of the complete global distribution of the z ~ 0 absorbers becomes necessary to 
resolve these issues. We have performed a comprehensive analysis of all the available data in the 
XMM-Newton archive. Our targets include a total of 43 AGNs. In this first paper, we present 
the results of our data analysis. We also report the first detection of z ~ 0 O VII absorption lines 
in 7 background AGNs at more than the 3cr level. Prior to this work, the most comprehensive 


study was performed by Miller Sz Bregman (2013). Their sample includes 26 AGNs, and 10 targets 
showed z ~ 0 O VII absorption lines at more than the 3<r level. Therefore our work is a significant 
increase in both the sample size and the newly detections of z ~ 0 O VII absorption lines. 


This paper is organized as follows. In section §2 we describe the observations and data analysis 
procedures. We present our main results in this section. We discuss the consistency of our data 
analysis by comparing our work with previous studies in section §3, as well as the impact of warm 
absorbers in some Seyfert 1 galaxies. We also comment on several targets with newly detected, 
z~0O VII lines in this section. The last section provides a summary. 


2. Observations and Data Analysis 


on-board 
on-board 
i collecting 


2.1. Target Selection 

Three high resolution X-ray spectrometers are suitable for the study of the narrow X-ray 
absorption lines: the low and high energy transmission gratings (LETG and HETG 
the Chandra X-ray Observatory, and the The Reflection Grating Spectrometer (RGS 
the XMM-Newton X-ray Telescope. We focus on the RGS data mainly due to its hig’ 
area at the wavelength region we are interested in. There are two RGS units: RGS 1 and RGS 2. 
However, due to operation failures the CCD 7 of RGS 1 and CCD 4 of RGS 2 are not working. 
Particularly, the CCD 4 of RGS 2 covers the He-like oxygen region which is our primary focus in 
this paper. Therefore in this work we will focus on the RGS 1 unit only. 

The z = 0 X-ray absorption lines from various ion species have been reported previously. These 
lines were produced by metal species in a variety of ionization stages, ranging from neutral to the 
H-like. Most reported ion species are Ne IX at 13.44 A, O VIII at 18.98 A, and O VII at 21.6 A 


(see, e.g, Nicastro et al.|2002 Fang et al.|2002 Rasmussen et al.|2003 Williams et al.|2005 ). In this 
work, we will mainly focus on highly ionized metals, in particular the He-like oxygen, O VII, for 
several reasons. First, Ne IX K a is located in a region where the RGS 1 has no effective area, and 
also there are several bad pixels near the O VIII I\ a region. Second, Under collisional ionization, 


1 See http://asc.harvard.edu/ 
2 See http://xmm.esa.int/ 
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the temperature for the peak O VII fraction is between 10 5 5 -10 6 " 5 degrees, providing an effective 
way to probe the hot gas in and around the Milky Way. 


We searched the entire XMM-Newton archive of the AGNs with enough RGS photon counts 
to warrant a spectral analysis. The RGS resolution has a full width of half maximum (FWHM) of 
~ 50 mA. We define one resolution element as half of the FWHM, or 25 mA. Our past experience 
indicates that the minimum requirement for a useful spectral fit is that the counts per resolution 
element, or CPRE, must be at least 20 photons. This selection criterion results in a total of 43 
AGNs in our sample. Among them 12 are BL Lac-type target. The remaining 31 are various 
Seyfert galaxies with a subclass of between 1 and 1.5. For convenience we broadly categorize them 
as Seyfert 1 AGN^J In Figure [l] we plot the all-sky Hammer-Aitoff projection of our targets in 
the Galactic coordinate. In Table [I] we list basic properties of these targets. Column (1) is the 
target name. Column (2) lists the AGN types. We broadly divide the entire sample into BL Lac 
targets (Type 1) and Seyfert 1 Galaxies (Type 2). Columns (3), (4), and (5) are the Galactic 
latitude, longitude and redshift, respectively. Column (6) is the Galactic neutral hydrogen column 
density (Dickey Sz Lockman 1990). Columns (7) and (8) are the total exposure time (in units of 
ksec ) and the CPRE. Columns (9) - (15) are line parameters which will be described in detail 
later. Columns (9), (10), (11) are the O VII column density, the Doppler-6 parameter, and the 
velocity shift of the line center, respectively. Columns (12) (13), and columns (15), (16) are the line 
line equivalent width (EW) and the signal-to-noise ratio (S/N) of the detection, for two different 
evaluation methods, respectively. We do not list the line parameters (denoted as “...”) if the line is 
not detected (defined as the S/N falls below 1), except the 3cr upper limits of the EW. We also list 
the C-statistic and degree of freedom in column (14) (see next section). In the last column we list 
the references that first reported the detection of the target. All the errors are la unless otherwise 
mentioned. We also plot the redshift distribution of our sample targets in Figure [2j The majority 
of the targets are located at z < 0.1. 


2.2. Absorption Line Analysis 

We followed the data reduction process using the standard XMM-Newton data analysis soft¬ 
ware SAS, version 12.0. iQ Specifically, we made use the tool rgsproc to process the RGS data. We 
rejected events with flags of “BAD_SHAPE”, “ON_BADPIX”, “ON_WINDOW_BORDER”, and 
“BELOW.ACCEPTANCE”. We also did not keep the cool pixels (“keepcool = no”). 

The above procedure outputs three files for each observation: A spectrum, a response file, and 
a background file. The reduced spectra were then analyzed with the X-ray software package XSPEC 


3 See the NASA Extrgalactic Database at http://ned.ipac.caltech.edu/ 
4 See http://xmm.esa.int/sas/ 
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Fig. 3.— RGS spectra of Mkn 421, PKS 2155-304, 3C 273, Ark 564, Mkn 509, and NGC 4051. 
The location of the O VII K a line is labeled with a vertical line. The red line in each panel is the 
model. The structures at 21.82 A are instrumental features. 


version 12.7.l|^J We fitted the continuum spectrum with a power law plus the absorption by the 
neutral hydrogen in our Galaxy. Since we are not interested in the broadband spectral properties, 
we selected only the spectral regions that contain the relevant spectral lines. The He-like oxygen 
has two transitions, the I\ a line at 21.6019 A and the Kp line at 18.6288 A (Verner et al. 1996). 
For most AGNs the K a line is the only prominent feature, so we performed the fit of the O VII K a 
line with a model that was developed in Buote et al. (2009) and discussed also in Fang et al. (2010). 
We ignored the regions below 21 A, and above 22 A, as well as the region around an instrumental 
feature at ~ 21.8 A. There are three free parameters in this Voigt-broadened line absorber model: 
the column density of the absorbing ion, the Doppler-6 parameter, and the redshift. This model 
reflects directly the physical properties of the absorbers and allows the line parameters to be tied 


5 See http://heasarc.nasa.gov/xanadu/xspec/ 
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Fig. 4.— Same as the Figure [3j but for RGS spectra for PKS 0558-504, MCG-6-30-15, Mkn 766, 
NGC 3516, H 1426+428, and MR 2251-178. 


consistently between the different spectra. The absolute wavelength uncertainty for RGS is about 
100 km s -1 0 so we limited the line center between with a redshift range of [-0.0004, 0.0004]. 
This range is also consistent with previous measurements in several high quality spectra. We also 
limited the Doppler-6 parameter in the range of 20 to 300 km s . The lowest temperature that 
can still produce a significant fraction of O VII is around 5 x 10 5 K under collisional ionization, 
this corresponds to a thermal velocity of ~ 20 km s” 1 . The upper limit is adopted by assuming the 
absorbing gas shall not escape from our Galaxy. We performed the fit by minimizing the C-statistic 
(identical to maximizing the Poisson likelihood function) of Cash (1979), which yields less biased 


best-fitting parameters (Humphrey et al. 2009). 


Figure [3 10 show the spectra of all our targets. Here again we emphasize that the spectra 


See XMM-Newton Users Handbook at http://xmm.esac.esa.int/external/xmm_user_support/documentation/uhb/index.html 
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Fig. 5.— Same as the Figure [3j but for Mkn 
and Mkn 279. 
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, Mkn 501, ESO 141-55, NGC 7469, NGC 3783, 


for each exposure were combined only for plotting purposes. The red lines are the best-fit model. 
In fitting the spectra of a target with multiple exposures, we did not co-add the spectra to form 
a single spectrum and response. Rather, we simultaneously fitted all the exposures, each with its 
own response. We tied together the three line parameters (column density, Doppler-6 parameter, 
and line center shift) among all the observation for the same target while we let the continuum vary 
for each observation. This procedure takes longer time to find the best-fit parameters, however, it 
avoids introducing systematic errors that may arise from co-adding spectra with different responses 
(see Rasmussen et al. 2007 for details). 


Among all the 43 AGN spectra in our sample, the O VII K a lines of 29 AGNs were discussed 
previously (see the last column of Table [I]). We have included 14 additional AGNs that satisfy 
our selection criterion. Table [l] columns 9-11 list the three fitting parameters and their la errors. 
Column (12) is the line equivalent width (EW), and (13) is the significance level. In total 21 targets 
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Fig. 6.— Same as the Figure [3j but for NGC 5548, H 2356-309, PKS 0548-322, NGC 4593, 
1H 0717-495, and Ark 120. 


show the O VII K a line at more than the 3a level, among which 7 were first detected in this worfQ 
We also list the C-statistic and the degree of freedom in column (14). In Table ?? we listed the 
detection of O VII Ka line at 3, 2, and la for different populations. 

For most AGNs we cannot constrain both the Doppler-6 parameter and velocity shift of the 
line center, so we also test another model in which we fix both values: We set the velocity shift of 
the line center at 0 km s _1 and Doppler-b parameter at 100 km s -1 . We listed the EW and the 
significance evaluated with this method in columns (15) and (16), respectively, in Table[l] While we 
see that both EW and the significance are slightly lower in this method, the results are essentiality 
consistent with those evaluated by allowing the line parameters to vary within a range. Unless all 


7 We emphasize that the RGS spectra of some of these targets have been analyzed before, but this is the first time 
that detections of the local O VII K a lines have been reported. 
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Fig. 7.— Same as the Figure [3j but for 1ES 1028+511, PG1116+215, 3C390.3, PG 1211+143, 
1H 1219+301, and PKS 2005-489. 


targets in actuality are properly described by these parameters, the derived line significances must 
be less in general, which is what we find. However, the differences are small and consistent with 
expectations that a range of 6-parameters and velocity shifts are present in the sample. 

Our default continuum model is a power law plus the Galactic neutral hydrogen absorption 
over approximately a 1 A range centered around the absorption lines of interest. The measured 
line properties crucially depend on the local continuum model, as well as the wavelength range 
over which the fitting is performed. To assess whether the power law model would be sufficient to 
characterize the continuum, we also investigated a more complex model, in particular a cubic-spline, 
to better match any fine details in the continuum. We found for most targets, such systematic errors 
are completely dominated by statistical errors. It only becomes important for the first several 
targets where the statistical errors are much smaller. In the case of Mkn 421, we found with a 
cubic-spline model the measured O VII I(a line EW is 12.47 ± 0.53 mA. This is slightly lower 
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Fig. 8.— Same as the Figure [3j but for Fairall 9, NGC 7213, Mkn 841, 3C 120, Mkn 110, and 
ESO 198-24. 


than what we measured using the single power law model. Such variations are consistent with the 
large-scale uncertainties in the RGS effective area on lA scale (Kaastr a et al.||20 06). We also tested 
the systematic errors that may be caused by the wavelength range of our fitting. We found that as 
long as we restricted our fitting to within a wavelength size of ~ lA, the uncertainties caused by 
systematic errors are negligible. 


A potential bias may occur when sometimes the fitting procedure tries to maximum the C- 
statistics. We tested our procedures by fitting simulated spectra. We define the bias as the difference 
between the input line EW and mean EW of the simulated spectra. We found in general for strong 
continuum with highly significant lines, the fitting procedure produced little or no statistically 
significant bias. For example, in the case of Mkn 421, the ratio between the bias and the standard 
deviation of the simulated EWs is less than 2%. When the continuum is weaker and the line is 
less significant, this ratio becomes slightly larger, but in general less than ~ 15%. So our fitting 
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Fig. 9.— Same as the Figure [3} but for IRAS 13349+2438, He 1143-1810, IC 4329a, TONS 180, 
Mkn 1502, and 3C 279. 


procedure shall not produce statistically significant bias. 


While our experience indicates that the such a model is often a reasonably good fit to the 
BL Lac-type targets, in the case of the Seyfert 1 galaxies, absorption features from the warm 
absorbers can sometimes substantially complicate the analysis, and each case has to be treated 
differently. We use NGC 3516 as an example here. In Figure [XT] we show the combined spectrum 
of NGC 3516 between 21.2 and 22 A. This target was observed six times between 2001 and 2006, 


with a total exposure time of 285 ksec (Mehdipour et al. 2010). Here in the plot the data and fit 


are for illustration purposes only since the actual fit was performed on the individual observation 
simultaneously. In the left figure we see the broad O VII absorption line produced by the warm 
absorber intrinsic to the AGN, centered at ~ 21.7 A. Using the power law plus the broad absorption 
line as our continuum model, In the bottom panel of the left figure, the Ay 2 plot clearly reveal 
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Fig. 10.— Same as the Figure [3j but for 1H 0414+009. 
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Fig. 11.— Spectra of NGC 3516 between 21.2 and 22 A before including a z ~ 0 O VII I\ Q 
absorption line (left figure) and after (right figure). The top panel in each figure is the spectrum, 
and the bottom is the Ay; 2 - The structure at rs_/ 21.81 A is an instrumental feature. The broad 
absorption feature at ~21.7 A is produced by the warm absorber intrinsic to NGC 3516. 


the narrow, z~0O VII absorption feature at around 21.6 A[^J In the left figure we show the fit 
after including the local O VII absorption line. The Ay 2 is improved substantially (A% 2 = 11 for 
77 degrees of freedom). 


Figure 12 shows the histogram distribution of the O VII I\ a Ka line EW (left panel) and 
column density (right panel). It appears that the most lines have an EW of ~ 20 mA, with a range 


of 10 - 30 mA. The column density distribution also centers around 
10 15 - 5 - 10 1 


10 16 cm 2 , with a range of 


,-2 


' 16,5 cm 2 . 


In Figure [I3| we show the dependence of the O VII absorption EW on the Galactic longitude 
(left figures) and latitude (right figures). Red squares and blue circles represent BL Lac and Seyfert 


8 While we use the C-statistics in this paper for less biased parameter fitting, here the conventional \ 2 statistics 
is adopted to demonstrate the significance of the absorption feature. 
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Fig. 12.— Histogram distribution of the O VII K a line EW (left panel) and column density (right 
panel). 


1 galaxies, respectively. For targets with a detection significance of at least la, we only plot the 
upper limits. To get a better understanding of the spatial distribution of the O VII absorbers, 
in the right panels of Figure [13] we divided the targets into two sectors based on their Galactic 


longitudes: Those close to the Galactic center within 45° (the top right panel of Figure 13) and 
those more than 45° away from the Galactic center (the middle right panel). The bottom left panel 
shows the dependence on the Galactic latitide of the entire sample. In the left panel of Figure [l3| 
we also divided the targets into two sectors based on their Galactic latitudes: Those close to the 
Galactic plane with |6| < 45° (the top left panel of Figure [l3|) and those that have \b\ > 45° (the 
middle left panel). The bottom left panel shows the dependence on the Galactic longitude of the 
entire sample. 

We have calculated the correlation between the line EW and the Galactic latitude and longi- 


tude, using software package ASURV Rev 1.2 ( 

Isobe & Feigelson 1990, 

LaValley et al. 1992 

), which 

implements the methods presented in 

Isobe et al. 

(1986 

). ASURV package is particular useful tool 


for correlation analysis of censored astronomical data, i.g., non-detections or detection limits (see, 
e.g., Keane et al. 2014). We listed in Table [3] the probability of correlation by chance for three 
correlation tests: Cox hazard model, generalized Kendall’s tau, and Spearman’s rho. A correlation 
exists if the probability is less than 5%. We have found in general the line EW has no correlation 
with either the Galactic latitude and the longitude. The only exception is the correlation between 
the line EW and the Galactic longitude for all the \b\ < 45° targets, in which all three tests give 
probabilities of less than 5%; however, a close look reveals that the correlation is largely dominated 
by one target, PKS 2005-398. After removal of this target, the probabilities increase to above 5%, 
and the correlation disappears. 


Gupta et al.| (2002| suggested that most lines they detected are saturated, leading to larger 
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Fig. 13.— 0 VII K a line EW as a function of the Galactic longitude (left panels) and latitude 
(right panels). Red squares and blue circles represent BL Lac and Seyfert 1 galaxies, respectively. 
For the Galactic latitude plot, the x axis is absolute latitude. For the Galactic longitude plot, the 
x axis is the angle from the Galactic center. The left panels from top to bottoms are for targets 
with \b\ < 45°, |6| > 45°, and all samples, respectively. The right panels from the top to bottoms 
are for targets within 45° from the Galactic center, more than 45° from the Galactic center, and 
all samples, respectively. 

ion column densities in their sample. The Voigt-fitting procedure we adopted allows us to extract 
directly the column density and Doppler-6 parameter information. Our study indicates that at least 
for the three brightest targets for which we can constrain the Doppler-6 parameter, the moderate 
variation of the Doppler-b parameters demonstrated by Mkn 421 (6 71 km s" 1 ), PKS 2155-304 

(~ 74 km s -1 ) and 3C 273 (~ 129 km s _1 ) does suggest a likely complicate environment for the 
OVII production. While opacity correction may be important for Mkn 421 and PKS 2155-304 due 
to their small 6 values, such correction shall become less important for targets such as 3C 273 with 
relatively larger 6-values. The line saturation may need to be quantified on a case-by-case base 
when allowed by photon statistics. 
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In Figure [14] we show the detection rate, or the sky covering fraction, as a function of the 
source flux. The sky covering fraction is defined as the ratio between the number of detections and 
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Fig. 14.— Detection rate as a function of the CPRE. The dark solid, red dotted, and blue dashed 
lines represent the detection rates for lines with at least 1, 2, and 3cr significance. The detection 
rates increase steadily for all the types to 100% for the brightest targets. 

the sample size. The solid black, dotted red, and dashed blue lines are for the detections with at 
least 3, 2, la significance, respectively. For the 3 a case, the covering fraction increases steadily 
from ~ 40% for targets with a CPRE of more than 20 counts, to 100% for targets with more than 
600 CPRE. For lcr significance, the covering fraction becomes 100% at ~ 100 CPRE. Our study 
suggests that the O VII absorbers have a rather uniform distribution in the CGM/ISM. 


3. Discussion 


3.1. Consistency Check with Previous Work 

The detection (or non-detection) of the O VII K a absorbers in a large fraction of our targets 
have been reported previously. Here we compare our results with those obtained previously to 
assess the consistency of our work. 

We first compare our results of the three brightest targets (Mkn 421, PKS 2155-304, and 
3C 273) with previous work since they have been studied multiple times by different groups, with 
both Chandra and XMM-Newton. We focus on the line EW since this is the observable that can be 
compared directly. We list the values of our work and several previous studies in Table |4j Columns 
2, 3, and 4 are EWs measured in this work, with XMM-Newton, and with Chandra , respectively. 
Our results are largely consistent with previous work within the error ranges. However, a noticeable 


difference is that for Mkn 421 and PKS 2155-304, our results are systematically higher than Williams 


et al. (2005) and William s et al.| (2007) by ~ 30%. While they used Chandra, for high quality data 
like these BL Lacs, such a discrepancy cannot be explained by the systematics between different 
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, in which 
. We refer 

the readers to this paper for further details. 


instruments ( Kaastra et al.|2006 ). Our results are consistent with Rasmussen et al. (2007 
they explained in detail the likely cause of this discrepancy with Williams et al. (2005) 


For the rest of the targets, we mostly compare with results from Bregman & Lloyd-Davies 


(2007), whose targets consist a subsample of ours. For most detections with more than la sig¬ 
nificance, our results are consistent with theirs within ~ 10% or within the error ranges. A few 
noticeable difference are as follows: In NGC 4051, PKS 0558-504, MR 2251-178, our best-fit values 
are different from theirs by a factor of 2-3; in NGC 7469 we have a 2.9cr detection while they had 
a non-detection; in PG 1116+215 we have a non-detection while they detected the O VII K a line 


at 2.2cr. Except for NGC 7469, our exposure times are significantly larger than those of Bregman 


& Lloyd-Davies (2007) (3 - 7 times longer), so our work should have a better statistical precision. 


For NGC 7469, we do not yet have a clear reason for the discrepancy. 


3.2. Galactic Absorbers or AGN Outflows? 


Highly ionized, X-ray absorbing gas has been frequently detected in the nearby Seyfert galaxies 
(see, e.g., Reynolds et al. 1997; Tombesi et al. 2013). While these systems offer an excellent 
opportunity to probe the inner environment of the AGNs, their presence poses a serious challenge 
to the study of the CGM/ISM in our Galaxy. These absorbers, in the form of warm absorbers 
or ultra-fast outflows, often have outflow velocities ranging from a few hundred km s -1 to mildly 
relativistic. Sometimes such a velocity can exactly cancel the redshift effect of the AGN, making 
the identification of these absorbers particularly difficult. A careful case-by-case study is therefore 
necessary to determine whether these absorbers are intrinsic to the AGNs or local in the CGM 
medium in our Galaxy. 


In our sample 31 background sources are Seyfert 1 galaxies, among which 20 show the detection 
of the local O VII K a line at more than the 2a level. We made a dedicated study of the warm 
absorbers reported in these 20 targets. We found 6 Seyfert 1 galaxies have no reported warm ab¬ 
sorbers; for the remaining 14 targets, we list them and information regarding their warm absorbers 
in Table [5] (also see the next subsection for comments on the targets with newly discovered z ~ 0 
O VII absorption line). Column 1 of Table [5] is the target name; column 2 is the recession velocity; 
column 3 is the reported velocities for various warm absorber outflows; and the last column is the 
reference in which the warm absorbers were discussed. We found that in NGC 4051 we do have 
difficulty separating the local absorber from various components of the warm absorbers. In the case 
of PG 1211+143, while Pounds et al. (2003) claimed that the O VII K a line is actually produced by 
a relativistic outflow with a velocity of ~ 24, 000 km s -1 , Kaspi & Behar (2006) found the absorbers 


can be explained as a much slower outflow but from other ion species. For all the other cases, it 
appears outflows cannot compensate the recession velocity, indicating the detected O VII K a lines 
are most likely local at z ~ 0. 
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There are further evidence that suggest most of the O VII absorbers should be local to our 
Galaxy. In a typical complex environment of the AGNs, the warm absorber/ultra-fast outflows are 
photo-ionized by the AGN, and exhibits absorption features in a variety of ion species in different 
ionization states. Sometimes emission lines are also detected. These absorbers/emitter typically 
coexist spatially, and form outflow components with different velocities (see NGC 7469 in our 
sample). However, for most reported O VII absorbers in our Seyfert 1 galaxy sample, we searched 
the entire spectrum for each data set and did not find other ion species (except in some cases Ne IX 
and O VIII) share similar outflowing velocity if assuming they are intrinsic to the AGNs. Such a 
single-line absorber is unlikely to be produced in the rich environment of an AGN. 


3.3. Comments on Individual Target with Newly Discovered Absorption Line 


We briefly discuss the targets that have a newly detected, z 
the 2 a level. 


0 O VII absorber at more than 


Mkn 335 — This is a Seyfert 1 galaxy at z = 0.0258. Previously neither X-ray or UV observa¬ 
tions showed any intrinsic absorption (see, e.g., Reynolds et al. 1997; Zheng et al. 1995). However, In 
2007 and later on, Mkn 335 went into a low state, providing an opportunity to study potential spec¬ 
tral features in its broad-line region (see, e.g, Grupe et al. 2007; Longinotti et al. 2008). Since then 
XMM-Newton observations revealed various outflow components at a velocity of ~ 5000 km s -1 
(see, e.g., Longinotti et al. 2013), while Hubble Space Telescope observations also indicated outflows 
in ultraviolet with veolcity of ~ 5000 km s" 1 (see, e.g., Longinotti et al. 2013). 


NGC 7 469 — This is also a Seyfert 1 galaxy at z = 0.0169 (de Vaucouleurs et al. 1991). This 


target was extensively studied for its ionized outflow and received broadband coverage in optical, 
ultraviolet and X-ray. Outflows at two main velocity regimes, 600 and 2300 km s -1 , were detected 


in both UV and X-ray observations (see, e.g., Scott et al. 2005 Blustin et al. 2007). Our detected 


feature at 21.6 A can also be interpreted as an outflow intrinsic to NGC 7469; however, that would 
imply a velocity of ~ 5000 km s -1 , which would be inconsistent with any of the known outflows. 
Furthermore, outflows often are associated with multiple features in emission and/or absorption. 
We examined the spectrum and cannot find any other lines from major ion species with similar 
velocity. Therefore, we conclude the feature at 21.6 A is most likely produced by O VII in our 
Galaxy. 


1H 0707-495 — This is a narrow-line Seyfert 1 galaxy at z = 0.0411 (Remillard et al. 1986). 


it is a well-studied target in the X-ray because of its broad iron emission lines that have been 
interpreted as relativistic emission lines from an ionized accretion disk (Fabian et al.| |2009). Lately 
a mild relativistic highly ionized outflow was also suggested (Dauser et al. 2012). The 21.6 A 


absorption feature was first noticed in Blustin & Fabian (2009); however, since they plotted the 
spectrum in the rest-frame of the AGN, this feature appeared at a wavelength of ~ 20.7 A, leading 
them to conclude this is an unknown feature. 
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PG 1211+143 - This is a bright quasar at z = 0.0809 (Marziani et al. 1996). Based on 
the first ~ 60 ksec observation with XMM-Newton, Pounds et al. (2003) reported the detection 


of a high-velocity ionized outflow with velocity around 24,000 km s i . Interestingly, as pointed 


out in McKernan et al. (2004), this velocity can cancel the redshift of the quasar and put the 


absorber in our Galaxy. Kaspi & Behar (2006) reanalyzed this observation and concluded that 
most absorption can be explained by a mild outflow with velocity around 3,000 km s^ 1 . This 
target was lately observed three more times with XMM-Newton. Here our analysis is based on 
all the four observations with a total exposure time of ~ 218 ksec. We confirmed the detection 
of the 21.6 A feature. We notice that in the Kaspi & Behar model, while they can identify most 


high-velocity ion species in Pounds et al. (2003) with a different ion species at a much lower velocity, 


they still cannot explain the features at 13.45 and 21.6 A. Coincidently, these are the rest-frame 
wavelength for the Ne IX and O VII K a lines, respectively, suggesting these two features may indeed 
be produced by the hot gas local to the Milky Way. 

ESO 198-24 — This is one of the two Seyfert 1 galaxies (z = 0.0453) in our sample that have no 
reported intrinsic absorption or warm absorber (Porquet et al. 2004). The only significant spectral 
feature that was detected was the Iron K a line at ~ 6.4 kev (see, e.g, Guainazzi et al. 2003). 

Mkn 84 I — This is also a Seyfert 1 galaxy at (z = 0.0364) that has been extensively studied by 
many X-ray satellites. Its soft excess can be explained either by the reflection from the photoionized 


disk (Crummy et al. 2006), or as ionized absorption in a relativistic wind (see, e.g., Gierlinski &; 


Done 2004 Petrucci et al. 2007). Longinotti et al. (2010) analyzed the high resolution RGS data, 
and found multiple component warm absorbers with velocities up to ~ 1600 km s -1 . 

Hell43-1810 — This is another Seyfert 1 galaxy (z = 0.0329) that has no reported intrinsic 


absorption or warm absorber. Cardaci et al. (2011) found a broad emission feature centering around 


21.45 A, with a width of ~ 2 A. They identified this feature as blended O VII emission triplets 
intrinsic to the source. 

1ES 1028+511, PKS 2005-489, and 3C 279 — All three targets are BL Lac-type objects 
in our sample. For 1ES 1028+511, previous Chandra observations with HETGS did not reveal 
any absorption/emission features (Fang et al. 2005); however, Steenbrugge et al. (2006) later re¬ 
ported tentative detections of narrow absorption features by intervening intergalactic medium, 
using Chandra LETGS and XMM-Newton RGS. No emission/absorption features were reported 
for PKS 2005-489 or 3C 279. 


4. Summary 

The z = 0 O VII absorption lines detected in the X-ray spectra of background AGN have 
a significant impact on our understanding of the Galactic gas content, both in the distant CGM 
environment or in the nearby Galactic disk. In this paper, we present an XMM-Newton survey of 
these local O VII absorbers. We summarize our findings here. 
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• We studied 43 AGNs that satisfy our selection criteria, 12 targets are BL Lac-type, and 31 
are Seyfert 1 Galaxies. A total of 21 targets show the detection of the z = 0O VII absorption 
line with more than the 3 a significance, among which 7 were newly discovered in this work. 

• We have fitted the z = 0 O VII K a transition with a Voigt-profile based line model. We find 
that most K a lines have an EW of ~ 20 mA, with a range of 10 - 30 mA. The column density 
distribution also centers around 10 16 cm 2 , with a range of 10 15,5 - 10 16 ' 5 cm 2 . 

• We correlate the line equivalent width with the Galactic coordinates and do not find any 
strong correlations between the line EW s with either the Galactic longitude or latitude. 

• The sky covering fraction, defined as the ratio between the number of detections and the 
sample size, increase from at about 40% for all targets to 100% for the brightest targets, 
suggesting a uniform distribution of the O VII absorbers. 

• Some AGNs have warm absorbers that may complicate the analysis of the local X-ray ab¬ 
sorbers. A case-by-case study indicates that while in very few targets we do have difficulty 
separating local and intrinsic absorbers, for the majority the 21.6 A absorbers appear to be 
located in our Galaxy. 

The study of the z = 0 O VII absorbers plays an important role in understanding the hot 
gas in and around our Galaxy, since O VII samples a broad range of temperature around 10 6 K. 
When combined with a study of the x-ray emission, an X-ray absorption study of the Galactic 
X-ray binaries, as well as a number of other investigations, we expect our study can help provide 
important clues to the properties of the CGM medium, the multi-phase model of the ISM, as well 
as questions such as the “missing Galactic baryons”. We plan to present our theoretical modeling 
in a later paper. 
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Table 1. The AGN Sample 
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Table 2. Target Detection 



> 3cr 

> 2(7 

> la 

Total 

This work 

7 

9 

11 

11 

BL Lac 

2 

3 

4 

4 

Seyfert 1 

5 

6 

7 

7 

All 

21 

28 

33 

43 

BL Lac 

7 

8 

10 

12 

Seyfert 1 

14 

20 

23 

31 


Note. — Statistics of the detected tar¬ 
gets. Top three rows show the detected tar¬ 
gets in this work; and bottom three rows 
show all the targets in our sample. Columns 
2, 3, 4 list the targets detected with a S/N 
larger than 3, 2, and ler, respectively. The 
last column is the total number of targets in 
each group. 



Table 3. Correlation Tests 0, 


Correlation 

Subsample 

Cox Hazard 

Kendall Tau 

Spearman Rho 

EW vs. Degree from GC 

6| < 45°, BL Lac 

6% 

7% 

10% 


6| < 45°, Seyfert 1 

14% 

29% 

22% 


6| < 45°, Total 

2% 

4% 

4% 


6| < 45°, Total b 

5% 

12% 

10% 

EW vs. Degree from GC 

6| > 45°, BL Lac 

86% 

90% 

88% 


6| > 45°, Seyfert 1 

93% 

89% 

94% 


6| > 45°, Total 

89% 

98% 

90% 

EW vs. Degree from GC 

BL Lac 

33% 

22% 

22% 


Seyfert 1 

44% 

47% 

46% 


Total 

25% 

18% 

22% 

EW vs. \b\ 

l < 45° or l > 315°, BL Lac 

11% 

12% 

16% 


l < 45° or l > 315°, Seyfert 1 

28% 

40% 

27% 


l < 45° or l > 315°, Total 

7% 

17% 

13% 

EW vs. \b\ 

315 > l > 45°, BL Lac 

46% 

83% 

87% 


315 > l > 45°, Seyfert 1 

49% 

96% 

94% 


315 > l > 45°, Total 

86% 

97% 

89% 

EW vs. \b\ 

BL Lac 

82% 

39% 

61% 


Seyfert 1 

16% 

56% 

56% 


Total 

24% 

33% 

42% 


Note. — a. Correlation tests using (1) Cox Hazard, (2) Kendall Tau, and (3) Spearman Rho methods. 
Columns 3—5 gives the probability that the correlation between eh two parameters is obtained by chance, 
and a less than 5% probability suggests the two parameters are correlated, b. Total sample excluding the 
bright BL Lac object PKS 2005-398. 
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Table 4. Comparison with Previous Observations 


Name 

This Work 

XMM-Newton 

Chandra 


(mA) 

(mA) 

(mA) 

Mkn 421 

13.13 ±0.67 

13.1 ± 1.0 [1] 

11-13 [2] 



11.8 ± 0.8 [3] 

9.1 ± 1.1 [4] 


PKS 2155-304 14.82 ± 1.20 13.7 ± 1.9 [3] 13.3 ± 2.8 [5] 

16.3 ±3.3 [6] 11.6 ±1.6 [7] 


3C 273 25.12 ±2.11 24.6 ± 3.3 [3] 28.4^ 2 2 5 [8] 

26.3 ±4.5 [6] 


Note. [1]: Rasmussen et al. 

(2007); [2] Kaastraet al. 

(2006); 

[3] Bregman & Lloyd-Davies 

2007); [4] [Williams et al. (2005); [5] 

Hagihara et al. |2010); [6] Rasmussen et al. 

(2003); [7] Williams 

et al. 

(2007); [8] Fang et al. ( 

2003) 
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Table 5. Seyfert 1 Galaxies with Warm Absrobers 


Name 

cz a 

Warm Absorber 

Reference 


(km s' 1 ) 

(km s' 1 ) 



Ark 564 7472 200-2000 

Mkn 509 10358 ~ 1000 

NGC 4051 800 200-600, 2340, 4600 

9000, 0.12c b 

MCG-6-30-15 2610 0, 1900 

NGC 3516 2201 ~ 100, 1000 

MR 2251-178 19194 < 130, ~ 480,15600 

Mkn 335 7740 5000 

NGC 7469 4619 600, 2300 

Mkn 279 9150 0-1500 

NGC 3783 2409 450-750 

1360-2080 


Smith et al. (2008) 


Kaastra et al. 


( 2012 ) 


Pounds & King (2013) 


Holczer et al. (2 

010) 

Mehdipour et al 

. (2010) 

Reeves et al. (2C 

)13) 

Longinotti et al. 

(2013) 

Blustin et al. (2 

307) 

Ebrero et al. (2C 

)10 

Kaspi et al. (20C 

>2) 


NGC 4593 

2651 

200-600 

Steenbrugge et al. 

(2003) 

1H 0707-495 

11880 

0.11-0.18c 

Dauser et al. 

(2 

012 



PG 1211+143 

24450 

3000 c 

Kaspi Sz Behar 

(2006 

) 



24000 c 

Pounds et al. 

(I 

>001 

1) 


Mkn 841 

10920 

< 1600 

Longinotti et a 

• (2 

010 

) 


Note. — (a): Recession velocity based on redshift. (b): c is the speed of 


light, (c): This is one warm absorber, but Kaspi & Behar (2006) and Pounds 
et al. (2003) have different interpretations. 





















































































































